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We report on the result of quantum Monte Carlo simulation of quasi-one-dimensional electron 
systems at 1/4-filling, considering organic superconductors such as TMTSF- and TMTTF-salts. 
We focus on the effect of dimensionality (interchain coupling) on superconducting fluctuation. 
First we consider Hubbard model which includes only on-site repulsion U. We find that the 
increase of interchain coupling enhances superconducting susceptibility, although it deforms 
the nested Fermi surface and suppresses the spin susceptibility. Next we consider an extended 
Hubbard model which includes nearest-neighbor repulsion V. In this case we consider the com- 
petition between different symmetries of electron pairing, d x 2_ y 2- and d^^-wave symmetry, and 
find that the particle-particle interaction vertex for d a .2„ H 2-wave pairing is no longer attractive 
even in the region where nesting condition still holds. On the other hand, the interaction ver- 
tex for d xy -wa,ve pairing persists to be attractive. The obtained results for the Hubbard model 
show qualitative agreement with recent renormalization group study, and have the possibility 
of accounting for the recent experimental result of (TMTTF)2SbF6, which exhibits a wider 
superconducting phase than the previous studies. 

KEYWORDS: superconductivity, quasi-one-dimension, extended Hubbard model, quantum Monte Carlo, 
organic conductors, TMTSF, TMTTF 



1. Introduction 

Quasi-one-dimensional (Q1D) electron systems, such 
as organic conductors TMTSF- and TMTTF-salts 
(TM2X), have been extensively studied both theoreti- 
cally and experimentally after the discovery of supercon- 
ductivity in (TMTSF^PFg. 1 Various experiments have 
been carried out and the results are summarized as the 
so-called Jerome's phase diagram. 2 In this generic di- 
agram, the ground state changes as spin-Peierls, spin- 
density-wave (SDW) and superconductivity (SC) by ap- 
plying pressure. The SC phase is located next to the 
SDW phase; it suggests that the source of the attractive 
interaction between electrons is due to antiferromagnetic 
spin fluctuation. The fact that the SC transition tem- 
perature, T c , decreases as applying pressure (going away 
from SDW phase) supports this prediction. Many the- 
oretical studies have been performed from this point of 
view 3 ~ 7 and the obtained results give qualitative agree- 
ment with the experiments. 

A recent experimental result exhibits, however, a wider 
SC phase under high-pressure. 8 This suggests that the 
effect of pressure (or increase of interchain hopping) on 
SC phase is positive. From the theoretical point of view, 
the possibility of higher T c due to the increase of inter- 
chain hopping has been proposed from the early days on 
the basis of one-dimensional (ID) renormalization group 
theory. 9 Recently, modification of this approach makes it 
possible to treat strong quantum fluctuations in Q1D sys- 
tem, 10 and concludes that the 6^2 ^2 -wave pair- field sus- 
ceptibility shows clear enhancement against interchain 
hopping. 11 ' 12 It indicates that the decline of T c along 
pressure may not be the case. 

However, the previous studies are mainly based on the 
perturbative methods, especially the ones which focus 
on the spin fluctuation. Although it has been verified 
that the results of these approaches qualitatively agree 



with that of numerical methods in two-dimensional sys- 
tems, 13 it has not been confirmed yet whether these ap- 
proximations are good in ID or Q1D systems. The main 
obstacle is the strong quantum fluctuations in these sys- 
tems, i.e., spin fluctuation and the other types of fluc- 
tuations are competitive and they interfere with each 
other. One of the methods which can treat these fluc- 
tuations correctly is the Q1D renormalization group the- 
ory mentioned above. However, it has its origins in the 
ID renormalization and is applicable only to wcakly- 
couplcd chains. Thus, in this paper, we adopt the quan- 
tum Monte Carlo (QMC) method in order to treat not 
only the spin fluctuation but also all types of fluctua- 
tions at the same weight, which is rather difficult in per- 
turbative approaches such as RPA or FLEX approxima- 
tion. In addition, it can treat any strength of interchain 
coupling without approximation in contrast to the Q1D 
renormalization group theory. We investigate the effect 
of interchain hopping, or dimensionality, on the SC fluc- 
tuation in the Hubbard model by using the auxiliary field 
QMC, 14 ' 15 having in mind the SC phase of TM 2 X salts 
under pressure. 

Furthermore, we study the effect of long range 
Coulomb interaction V, which has been recognized to 
be very important in organic conductors, especially the 
TMTTF-salts, in the context of the charge ordering 
(CO). 16 ' 17 So far, the effect of V on the SC has been in- 
vestigated for the SC phase next to the CO phase, such 
as a-ET 2 I 3 , 18 - 20 or /?-(DMeET) 2 PF 6 . 21 As for TM 2 X 
salts, on the other hand, it has not been studied since 
the SC phase is separated from the CO phase. However, 
the high-energy peak (~ 30meV) of the optical conduc- 
tivity has been observed even in the metallic TMTSF 2 X 
salts, 22 indicating a precursor of CO. The effect of V, 
therefore, is not negligible even in the metallic QID com- 
pounds; it can give an sizable impact on the mechanism 
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of SC. In order to clarify this effect, we also investigate 
the behavior of SC susceptibility for the Q1D extended 
Hubbard model. 

This paper is organized as follows. In §2, we introduce 
the model and method used in this paper. In §3, the 
results of the Hubbard model are presented, which cor- 
respond to the situation similar to the previous studies. 
They show qualitative agreement with the renormaliza- 
tion group study. 12 In §4, the result of the extended Hub- 
bard model is presented. The effect of V mainly appears 
as the change of SC symmetry from d x 2_ y 2 to d xy . The 
comparison with the experimental results is given in §5. 

2. Model and Method 

We consider the two-dimensional Hubbard and ex- 
tended Hubbard model on an anisotropic square lattice, 
whose Hamiltonian is given by 

H = ^ {-Uj4 a c ja + h.c.) + U J] n iT n a 

(ij)" « 

+ v 53 n i" n 3"' - A f 53 nm ' w 

(ij)cr<T f i& 

where cj CT (cj CT ) is the creation (annihilation) operator 
of an electron with spin a (a =T,I) on site i and 
Uia = c ia Ci a . Here tij has two different values, t x (intra- 
chain hopping) and t y (interchain hopping) . The nearest- 
neighbor Coulomb interaction V is introduced as de- 
picted in Fig. 1. We set the intrachain hopping t x = 1, 



We calculate spin, charge and pair-field susceptibilities 
Xs(q) = ^ f dr^e^-^KCrKtO)), 

ij 

Xc (q) = 1 [ drJ^^-^HniWhM), 
iV Jo ij 

X d sc = ^ / /3 dr53(A J (r)At(0)) (2) 

ij 

with 

m i( T ) = 4^ T ) c n( T ) ~ c li( r ) c u( r ). 
fti(r) = c\^{t) Ci] {t) +4^)0^) - (n), 

A *( T ) = ^2fd($)Cit(T)Ci + 6l(T). 
S 

Here fd{5) is the factor +1 or —1 corresponding to the 
pairing symmetries as shown in Fig. 2. 



+ 




Fig. 2. The internal pair coordinate 8 and the factor fd{S) used 
in defining Ai for d x 2_ y 2-v/ave pairing (left panel) and d xy -w&vc 
pairing (right panel). 
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Fig. 1. Transfer integrals t x 
V. 



t y and nearest-neighbor repulsion 



interaction strength {7 = 4 and k^, = h = 1 through- 
out this paper. The number of electrons are set to be 
(n) = 1/4 by tuning the chemical potential fi. 

In order to investigate the properties of this model, 
we employ the auxiliary field quantum Monte Carlo 
method. 14, 15 The shape of the cluster is set to be 16 x 8 
for V = or 16 X 4 for V ^ 0. It means that there are 
16 sites along chains (x-direction) and 8 (or 4) chains 
are coupled with each other in the y-direction. In the 
presence of V, we restrict the system size to 16 x 4 and 
the temperature range to /?(= l/T) < 4, because of 
heavy computational cost of Monte Carlo sampling. 25 
The imaginary-time step is set to be At = 1/8, and 
64000 measurements separated by two sweeps are per- 
formed after 1000 warm-up. 



In addition to eq. (2), we also calculate the uncorre- 
cted pair-field susceptibility 26, 27 



XdSC 



J0 iSjS' 



x (c^r^mia+siir)^,^)). (3) 

-G,,i,tW -G i+S:j+S , u (r) 

The difference between eqs. (2) and (3) is whether the 
particle-particle vertex T is included or not (see Fig. 3). If 
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Fig. 3. Full susceptibility XSC & n d uncorrelated susceptibility 
Xsc- Thick lines denote the dressed Green's functions. T is the 
particle-particle interaction vertex. 



Xsc is larger than xsc, it means that T is attractive. 26-28 
We also calculate the density of states -ZV(w), which is 
defined as 



G n {r) = - 



dcu- 



-/3uj 



(4) 
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where GV,(t) is the imaginary-time Green's function 
Gti(T)=-£<Mr)4(0)>. 

This expression will be readily understood from the def- 
inition of N(u>) 

k 

with A{k,uj) being the one-particle excitation spectrum 



G(fe,r) = - 



dcu- 

1 4- c _/3w 



The inverse transformation from Gu(t) to N(u>) in eq. 
(4) is carried out by the maximum entropy method com- 
bined with QMC. 29 

3. Results of the Hubbard model 

First we show the results of the Hubbard model, whose 
Hamiltonian is given by eq. (1) with V = 0. Figure 4 is 
the spin susceptibility Xs{Q) an d SC susceptibility Xdsc 
against temperature T for different values of t y . As for 
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Fig. 4. Temperature dependences of spin susceptibility Xs(Q) 
(upper panel) and d a! 2_ !/ 2-wave pair-field susceptibility XdSC 
(lower panel) for U = 4 and different values of t y . The change 
of these susceptibilities when t y is increased are indicated by ar- 
rows. Where not shown, error bars are smaller than the symbols. 
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Fig. 5. Fermi surfaces with t y = 0.0 (solid), 0.2 (dotted) and 1.0 
(dash-dotted) at 1 /4-filling. The good nesting condition holds at 
about 0.0 < t y < 0.3. Q = (t/2, it) is a nesting vector. 



ID, but the enhancement becomes moderate as t y in- 
creases. This is due to the breaking of the nesting condi- 
tion, as depicted in Fig. 5. On the other hand, Xdsc shows 
different behaviors. Xdsc is enhanced by t y in small t y re- 
gion, while Xs{Q) decreases. XdSC becomes largest at the 
intermediate values of t y , such as t y = 0.4 or 0.6. 

The iy-dependences of Xs(Q) an d XdSC ar< 5 summa- 
rized in Fig. 6 when the temperature is fixed at T — 
0.167. The non-monotonic behavior of XdSC can be in- 
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Fig. 6. ty-dependences of spin susceptibility Xs{Q) (upper panel) 
Xs{Q), it is Strongly enhanced at low temperature near and ^2.^2 -wave pair-field susceptibility \dSC (lower panel, solid 

squares) at T = 0.167 for [7 = 4. The open squares denote the 
uncorrelated pair-field susceptibility XdSC- 
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tcrprcted as follows: the enhancement comes from the 
pair-hopping of electrons, and the suppression is due to 
the suppression of spin fluctuation because of the defor- 
mation of the nested Fermi surface. In fact, the peak posi- 
tion of Xs(q) locates at around q = Q for 0.0 < t y < 0.6, 
while it does not at t y = 0.8 and 1.0 (not shown). 

In Fig. 6, the uncorrelated susceptibility, XdSC, is also 
shown in comparison with the full susceptibility, XdSC- 
We find that the vertex part is attractive for 0.0 < t y < 
0.6, while it is repulsive for t y — 0.8 and 1.0. This is 
consistent with the above picture that XdSC is most en- 
hanced at the region where the nesting property and the 
interchain pair-hopping of electrons arc well balanced. 

In order to confirm the effect of nesting property, it is 
useful to introduce the frustrated hopping, t', as shown in 
Fig. 7. This kind of hopping is another key to control the 



, t' 



Fig. 7. Left panel: the interchain hopping t', which induces geo- 
metrical frustration. Right panel: Fermi surfaces with frustration 
t' when t y = 0.0 (solid), ty = 0.2 (dotted) and 1.0 (dash-dotted) 
at 1/4- filling. The value of t' is chosen as t' = t y /2. It should be 
noted that the nesting vector Q no longer exists. 



strength of spin fluctuation and interchain coupling. In 
fact, this type of transfer integral has been confirmed to 
exist in many TM^X-salts. 30 It deforms the Fermi surface 
and worsen the nesting condition (compare Fig. 5 and 7). 

Figure 8 shows the ^-dependence of spin and SC sus- 
ceptibilities at T = 0.167 for different strengths of t'. 
Here we choose the values of t' as t' = 0, t y /2 and t y . It 
is seen that Xdsc is slightly enhanced at small t y region 
and suppressed at large t y region, while Xs(Q) is mono- 
tonically suppressed by t' . The enhancement of Xdsc will 
be again due to the increase of interchain coupling. In- 
deed, the interchain hopping t' couples the chains more 
tightly, while it worsen the nesting condition. It will in- 
crease the pair-hopping of electrons and, as a result, XdSC 
is slightly enhanced by t' . In addition, the region where 
the interaction vertex is attractive (xdsc — Xdsc > 0) 
shifts to the smaller t y region (not shown), which simply 
supports this picture. 

However, one might point out that the maximum of 
Xdsc is due to the van Hove singularity which exists at 
t y ~ 0.35 in 1/4-filled non-interacting system, i.e., the 
effect of electron-electron interaction does not matter. 
In order to clarify this point, we calculate the density 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



t'=0 




0.2 0.4 0.6 0.: 




Fig. 8. t y -dependences of spin susceptibility Xs(Q) (upper panel) 
and d x 2 _ y 2-wave pair-field susceptibility XdSC (lower panel) at 
T = 0.167 for U = 4. Squares, circles and triangles correspond to 
the case tl =0, ty/2 and ty. The change of these susceptibilities 
when t' is strengthened are indicated by arrows. 



of states (DOS) N(ui) for interacting system. Figure 9 
shows N(ui) at T = 0.167 for different values of t y (here 
we consider the case t' = 0). Although the suppression 
of N(u> = Li) at large t y region is recognized, the en- 
hancement at small t y region seems not to exist. This 
behavior does not follow the curve of Xdsc against t y . 
Thus we conclude that the effect of van Hove singularity 
is irrelevant to understand the non-monotonic behavior 
of XdSC: a s long as the temperature range we calculated. 

4. Results of the extended Hubbard model 

In the previous section, we investigate the case with 
V = 0, where Xdsc is affected by the pair-hopping of 
electrons, in addition to the strength of spin fluctuation. 
On the other hand, the pairing mechanism itself is well 
understood by the spin-fluctuation theory. In this sec- 
tion, we consider whether this pairing mechanism is valid 
or not under the presence of charge fluctuation induced 
by the nearest-neighbor repulsion V. 

We consider the extended Hubbard model, whose 
Hamiltonian is given by eq. (1). We impose the limitation 
on not only the system size and temperature but also V 
as V < 0.5 in order to avoid round-off errors and severe 
sign problem. (We adopt the parameters which assure 
that the average sign (S) becomes greater than 0.1.) 

First we show the spin and charge susceptibilities with 
and without V. The upper two panels of Fig. 10 rep- 
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Fig. 9. Density of States for (7 = 4 and different values of t y at 
T = 0.167. The frustration t' is not included. 

resent the spin susceptibility, x s (q), for V = and 0.5 
when T = 0.25 and t y = 0.2, where the temperature is 
relatively low and the nesting condition still holds. It is 
seen that the structure of Xs(q) is hardly affected by V. 
It comes from the fact that Xs{q) is mainly affected by 
the shape of the Fermi surface, which is considered not 
to be drastically deformed by V. The lower two panels 
of Fig. 10 show the charge susceptibility, Xc(o)- A sharp 
peak at q = (tt, tt) is formed when V is finite, which is 
considered to be the reflection of the lattice structure. 
(Note that we consider the 1/4- filled system, which fa- 
vors q = (tt, tt) charge disproportionation under finite 
V".) This result suggests that the charge fluctuation in- 
duced by V is more affected by the geometrical configura- 
tion of V rather than the shape of the Fermi surface. 31 It 
is a clear contrast between Xs and Xc : the nesting prop- 
erty (momentum space) and the lattice structure (real 
space). 

It is an interesting subject to investigate how V affects 
SC susceptibility. To understand the effect of charge fluc- 
tuation, we consider not only the d x 2 _ y 2-wave pairing 
but also d xy -wa,ve pairing, whose gap function Aj is de- 
fined in Fig. 2. Figure 11 shows the temperature depen- 
dence of the SC susceptibilities for these two pairings. 
They are qualitatively similar, but it seems that d xy - 
wave pairing is slightly larger than the d x 2_ y 2-w&ve pair- 
ing. This advantage of d xy -wsive pairing over d x 2_ y 2--wave 
pairing will be recognized more clearly in Fig. 12, where 
the independences of these susceptibilities at T = 0.25 
are shown. (The results without V are also displayed 
for comparison.) The d x 2_ y 2-wsve pairing is drastically 
suppressed by V, although the suppression of d xy -w&ve 
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Fig. 10. Spin susceptibility Xs(<j) (upper panels) and charge sus- 
ceptibility Xc(<?) (lower panels) when T = 0.25 and t y = 0.2 
in the 1st Brillouin zone. The interaction strengths are U = 4, 

V = (left panels) and U = 4, V = 0.5 (right panels). The struc- 
ture of Xs(<?) reflects the shape of the nested Fermi surface with 
a nesting vector Q = (n/2,n) and seems not to be affected by 
V. On the other hand, the structure of Xc(<?) is largely affected 
by V and its peak position reflects the geometrical structure of 

V (or its Fourier transform V(q) oc cos(q x ) + cos(q y )). Note that 
each point includes errors with a few percent. 

pairing is not so large. Furthermore, it should be noted 
that, for V = 0.5, the interaction vertex for d a ,2_ J/ 2-wave 
pairing is repulsive for any value t y , even for t y = 0.0 or 
0.2 where strong spin fluctuation still exists. This sug- 
gests that the spin and charge fluctuations cancel each 
other and the vertex part for d x 2_ y 2-w&ve pairing is no 
longer attractive. On the other hand, the interaction ver- 
tex for the d X y-wave pairing is still attractive; it is even 
enhanced at the region of large t y . 

Finally, we estimate the effects of V on the dynamical 
property of this model. Figure 13 shows the DOS, N(oj), 
for different values of t y at T = 0.25 and V = 0.5. Note 
that N(lo = fi) on the Fermi surface is finite for any 
value of t y , i.e., the system is metallic and the CO does 
not occur. It is surprising that such a small V largely 
affects the SC interaction vertex. 

In this section we have investigated the effect of charge 
fluctuation induced by the nearest-neighbor repulsion V 
on the SC symmetry. Our results indicate that the sym- 
metry of electron pair changes from d x 2_ y 2 to d xy , even if 
there still exists strong spin fluctuation and V is as weak 
as not accompanying CO. It suggests that the effect of 
V on SC is not negligible, even in the metallic system 
with strong spin fluctuation. We should incorporate not 
only the shape of the Fermi surface but also the effect of 
V(q) into discussion. 
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Fig. 11. Temperature dependences of SC susceptibility for 
_ !/ 2-wave (upper panel) and d xy -wave (lower panel) pairings 
for (n) = 0.5, L7 = 4 and V = 0.5. 



5. Conclusions 

In this paper we have investigated the superconduc- 
tivity in Q1D system from the non-perturbative point of 
view, focusing on the effect of dimensionality on SC fluc- 
tuation. First we introduced the result of the Hubbard 
model, which includes only the on-site repulsion U. In 
this case, it is found that the d a .2_, y 2-wave pair-field sus- 
ceptibility, Xdsc, is increased by t y in the small t y region, 
while the spin fluctuation Xs(Q) is monotonically sup- 
pressed by t y . This behavior comes from the fact that the 
increase of t y has two different effects; the deformation of 
the Fermi surface and the enhancement of pair-hopping 
of electrons. The d x i _ y 2-wave SC is most enhanced at 
the region where these two factors are well balanced, as 
shown in the left panel of Fig. 14 schematically. This 
behavior suggests that the monotonic decline of T c may 
not be the case. On the other hand, the pairing mech- 
anism itself is the well-known spin fluctuation assumed 
by many previous studies. It is confirmed by the fact 
that the shape of the Fermi surface determines whether 
the particle-particle interaction vertex T for d x 2_ y 2-wave 
pairing is attractive or not. 

Next we considered the effect of nearest-neighbor re- 
pulsion V, focusing on the effect of charge fluctuation 
induced by V on the symmetry of electron pair. The ex- 
pected behavior for d xy -wave pairing is shown in the right 
panel of Fig. 14. We speculate that the charge fluctuation 
induced by V is not affected by the nesting property of 



Fig. 12. (Color online) t^-dependence of SC susceptibility at T = 
0.25 for U = 4, V = 0.0 (upper panel) and U = 4, V = 0.5 
(lower panel). The solid and open symbols denote the full and 
uncorrclated susceptibility. The contribution of particle-particle 
interaction vertex for d x 2_ y 2-wa,ve pairing is indicated by arrows: 
attractive (upper panel) and repulsive (lower panel). 
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Fig. 13. Density of states at T = 0.25 for U = 4 and V = 0.5. 
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Fig. 14. Schematic representation of the effect of dimensionality 
on SC fluctuation in Q1D 1/4-fillcd system. 



the Fermi surface but the geometrical lattice structure. 
In this case, the increase of t y simply makes a positive 
contribution to SC because of the enhancement of the 
electron pair-hopping. Our results support this picture, 
and in addition, indicate that T for _ y 2 -wave pairing 
is repulsive when V is finite, even at the region where 
strong spin fluctuation still exists and V is as small as not 
accompanying CO. On the other hand, d xy -wave pairing 
is not so largely affected by V and T for d^-wave pairing 
remains attractive. This indicates that the nesting prop- 
erty is not necessarily the crucial factor for determining 
the symmetry of electron pair, even if the system has a 
good nesting property. It may sound peculiar, especially 
from the perturbative point of view which focuses on the 
pairing mechanism mediated by the fluctuation of the 
ordered phase next to the SC. However, the possibility 
of strong cancellation between spin and charge fluctua- 
tions in r would not be so unrealistic. It suggests that 
we should incorporate the effect of long range Coulomb 
interaction V into the discussion of SC in order to deter- 
mine the symmetry of electron pair. 

Let us compare our results with experiment. We 
showed that the Xdsc is enhanced by increasing t y , 
whichever the symmetry of electron pair is d x i_ y i or d xy . 
It suggests that the SC phase in TM 2 X salts have the 
possibility of becoming wider, such as what is confirmed 
in (TMTTF) 2 SbF 6 . 8 Furthermore, we also found that 
the dry-wave pairing becomes more realistic than d x i_ y i- 
wave pairing under the presence of nearest-neighbor re- 
pulsion V. It indicates that the d x 2_ y 2-w&ve pairing 
which is assumed in the previous studies of TrV^X-salts 
may need to be revised. In fact, not only the TMTTF- 
salts, the TMTSF-salts showing the metallic behavior at 
ambient pressure have the precursor of CO as the higher- 
energy peak of the optical conductivity. 22 This means 
that the long-range Coulomb interaction is not negligi- 
ble in not only TMTTF-salts but also TMTSF-salts. 

Finally, let us comment on the temperature range of 
the present calculation. We cannot study very low tem- 
perature region due to the severe sign problem. This 
could make trouble with the analysis of interaction vertex 
T; it is possible that the strengths of the spin and charge 
fluctuations which contribute to T changes as tempera- 
ture is lowered, since the spin susceptibility Xs(Q) grows 
rapidly while the charge susceptibility Xc(tt, tt) shows 
more moderate behavior in the temperature range we 
calculated (not shown). Thus, the effect of V on SC sym- 



metry at lower temperatures still remains open. We hope 
that more studies would be carried out both from theo- 
retical and experimental point of view. 
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the shape of the Fermi surface. However, the effect of V(q) V which reflects the lattice structure rather than the nesting 

is stronger than that of Xo(<?)- This leads to the emergence condition, 
of q = (ir, it) charge disproportionation under the presence of 



